A new cooperation mechanism of kinesin motors when extracting membrane tube 
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Membrane tubes are important elements for living cells to organize many functions. Experiments 
have found that membrane tube can be extracted from giant lipid vesicles by a group of kinesin. 
How these motors cooperate in extracting the fluid-like membrane tube is still unclear. In this 
paper, we propose a new cooperation mechanism called two-track-dumbbell model, in which kinesin 
is regarded as a dumbbell with an end (tail domain) tightly bound onto the fluid-like membrane 
and the other end (head domain) stepping on or unbinding from the microtubule. Taking account 
of the elasticity of kinesin molecule and the exclude volume effect of both the head domain and the 
tail domain of kinesin, which are not considered in previous models, we simulate the growth process 
of the membrane tube pulled by kinesin motors. Our results indicate that motors along a single 
microtubule protofilament can generate enough force to extract membrane tubes from vesicles, and 
the average number of motorspuUing the tube is about 8 ~ 9. These results are quite different 
from previous studies (Ref. [131), and further experimental tests are necessary to elucidate the 
cooperation mechanism. 

PACS numbers: 87.16.Nn 



I. INTRODUCTION 

Membrane tubes widely exist in eukaryotic cells, such 
as in the endoplasmic reticulum and the Golgi apparatus 
etc.[l|, and play an important role in intracellular trans- 
portation as a long-range tubular transport intermediate 
which is essentially different from the conventional no- 
tion of small vesicle transportation [2]. Experiments have 
revealed that the membrane tube in living cells always 
co-occur with microtubules and recent experiments 
in vitro have shown that both the processive and the 
non-processive motors can extract membrane tubes from 
lipid giant unilamellar vesicles 0411 ■ Theories and exper- 
iments also showed, after the formation of a membrane 
tube from a vesicle, the force needed to keep the growth 
of the membrane tube is about /o ~ 27.5 ± 2.5 pN [6, 9], 
which is far greater than the stall force of a single motor, 
about 6~7 pN So the extraction of membrane 

tube must be a result of cooperation of many motors. 

Unlike the widely-studied rigid or elastic cargo 
membrane tube are actually fluid like. The lipids on the 
membrane tube, which the motors are adhered to, can 
flow on the surface of membrane tubes. So some models 
suggest that only the leading motor can apply force to 
the membrane tube [l^ [l^. According to these mod- 
els, however, it is obvious that the motors can not ap- 
ply enough force to extract membrane tube from vesi- 
cles. Campas et al. have considered various motor coop- 
eration schemes (see Ref.fT?] for details) and suggested 
that there are three protofllaments of the microtubule 
simultaneously involved in the pulling process. This pro- 
posal seems plausible, but no experiments support it di- 
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rectly. Especially, the structure and elasticity of motors 
are not considered in the former models, while these de- 
tails can substantially affect the collective transport of 
motors |14il . So the cooperation mechanism should be 
further revised, and other mechanism could be possible. 

In this paper, we focus on the extraction mechanism of 
membrane tube by processive motors (kinesin) and pro- 
pose a two-track-dumbbell model (see Sec. HI] and Fig. [T] 
for details ) , taking account of the elasticity of the motor 
molecule and the excluded volume effect of both the head 
domain and the tail domain of the motor molecule. By 
this model, we conclude that motors along a single mi- 
crotubule protofllament can also generate enough force 
to extract membrane tubes from giant vesicles. 

II. MODEL 

In the experiments of Ref. [l3l , each kinesin binds to 
a rhodamin-labeled biotinylated lipid on the vesicle mem- 
brane through a streptavidin molecule. And the density 
of kinesin on the membrane can be directly controlled by 
flxing the biotinylated lipid concentration in the mem- 
brane. It has been found that if the density of kinesin 
on the vesicle exceeds a threshold about 100~200 ^m^^, 
membrane tubes can always be formed; and if the den- 
sity is smaller than 100 /im~^, no tube extraction can be 
observed over a long time period. These experiments sug- 
gest that if the density is low, there are no enough motors 
gathering at the tip of membrane tube, and no enough 
force can be generated to maintain the tube growth. So 
a group of kinesin on the tip is essential for the formation 
of membrane tube. As mentioned in Sec. U however, it 
seems that only the leading one of the motors can ap- 
ply force to the membrane tube, but generate no enough 
force to resist the retreat of membrane. So there must 
exist some mechanisms, in which the non-leading motors 
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at the tip region of membrane tube can also apply force 
to the membrane tube. 

Noticing the finite volume of the complex of motor tail 
domain and streptavidin. The tail domain complexes of 
neighboring motors may exclude each other, and generate 
repulsive force in between. In this way, the non-leading 
motors can also apply force to the leading motor and thus 
to the membrane tube. On the other hand, the elasticity 
of the stalk of kinesin molecule has been found to play an 
interesting role in multi-motor transportation in our pre- 
vious study [ij]. In the present paper, the kinesin stalk 
is also considered as a spring, with an elastic constant 
k 



ward with a rate fc^ or detach from the microtubule with 
a rate 



spring • 



Combining the above two points, one can pic- 
ture the kinesin molecule as a dumbbell with a spring-like 
stalk (Fig. 1). One end of the dumbbell (the tail domain) 
can slide on the surface of the membrane tube, and can 
not detach from membrane tubes because of strong con- 
nection between kinesin and lipid. The other end of the 
dumbbell (the head domain) can either bind to or unbind 
from the microtubule, and the bound motor can step for- 
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FIG. 1: (Color online) The two-track-dumbbell model with 
consideration of the elasticity of motor stalk. The membrane 
tube is discretized as one-dimensional lattice, with the the 
same lattice spacing as that of the microtubule. The num- 
ber in each box represents the number of free motors in each 
membrane site, (b) shows the situation after the leading mo- 
tor unbinding from microtubule. Free motors are dragged 
back by the membrane tube and change their positions ac- 
cordingly upon the motion of membrane tube, (c) shows the 
situation after the leading motor stepping forward. Free mo- 
tors are dragged forward by the membrane tube. The number 
of free motors at site 1, mi, adopts 1 with a possibility of mo, 
else mi = 0, which flows from the vesicle. In (b) and (c) 
the actual length of the membrane tube is determined by the 
equilibrium of force, while the figure just to give a demonstra- 
tion. 



ku{f) = fc«oexp(///d) 



(1) 



where fc„o is the unbinding rate at vanishing load, and 
Jd is the detachment force with the value of about 3.0 
pN The backward stepping of motor is neglected 

here. The motors which attach to the membrane tube 
but not bind to the microtubule (i.e., free motors), can 
diffuse freely on the surface of membrane tube with a 
diffusion constant D and bind randomly at a rate fcf, to 
non-occupied sites on the microtubule. 

To carry out a simulation of the tube growth (Sec lIIII) . 
we have adopted a discrete microscopic approach in 
which the membrane tube is discretized into M + 1 sites 
(called membrane sites) numbered from to M (Fig. 
1), and the lattice spacing d equals to the periodicity 
of microtubule, 8nm. The boundary membrane site is 
connected to the reservoir of motors (the lipid vesicle) 
which can supply motors continuously when membrane 
tube is growing, and the motor density at this site is po, 
i.e. the motor density on the vesicle. So the number of 
motors on membrane site is mo = 2iTrdpo, where the 
radius of membrane tube r is related to the retreating 
force F and the bending modulus k of the membrane by 
r — 2ttk/F 9]. Since the whole perimeter of the tube 
exceeds 60 nm, each membrane site can contain many 
motors. In our simulations, however, we find that the 
number of free motors at a given membrane site seldom 
exceeds two, which means the free motor density is so 
low that their diffusion is almost uncorrelated. So the 
diffusion rate of a free motor along the membrane tube 
can be taken as kd — D/d^ 17]. A bound motor steps 
forward to the next un-occupied microtubule site with 
a rate kf{f) = V{f)/d, where V{f) is the force-velocity 
function for single motor transport which has been widely 
studied both experimentally and theoretically [13, Ell- In 
this study, the V{f) function is adapted from theoreti- 
cal results of two-state stochastic model of Fisher and 
Kolomeisky, with [ATP]=1.0 mM (see Ref. [lH for de- 
tails) . 

We now turn to the tip of membrane tube. Each bound 
motor behaves just like a spring connecting two beads. 
One bead (the head domain) binds to the microtubule, 
while the other (the tail domain) connects tightly to the 
membrane tube. If a motor bears force, the spring is 
stretched. In this model, the loading force is directly 
transmitted to the leading bound motor. The following 
motors can also share the load, by pushing the leading 
motor due to the exclusion of their tail domains (Fig. [!}. 
In fact, most of the tip motors would bear the load. For 
instance, head domains of two neighboring motors may 
be separated by one or more empty microtubule sites, 
but their tail domains can still interact with each other 
since their stalks can be stretched (see Fig. [TJa), the 
five leading motors on the tip will share the load). Mo- 
tors which bear the load are called pulling motors. The 
length of membrane tube is determined by the equilibra- 
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tion between retreating force and extracting force which 
depends on the distribution of bound motors on the mi- 
crotubule. 

If a motor unbinds from the microtubule, it can diffuse 
on the surface of the membrane tube (Fig. [Ijb) ) . If a mo- 
tor steps forward on the microtubule, its tail domain will 
be dragged forward (Fig. [lie)). For an unpulling motor, 
the tail domain is usually dragged forward by one mem- 
brane site after a forward step, and it does not change 
the tube length and the state of other motors. But an 
unpulling motor can also become a pulling motor when 
it encounters another pulling motor after a forward step. 
Then the load should be redistributed among pulling mo- 
tors. 



III. SIMULATION ALGORITHM 

We adopt the Gillespie algorithm [IS] to simulate the 
growth process of the membrane tube, based on the 
model introduced in Sec. [ill 

As mentioned, each bound motor can unbind from mi- 
crotubule and step forward, and the free motors can dif- 
fuse and bind to microtubule. In the simulation, we de- 
fine an array Pi numbered in order I ^ [1, 2N + 3M + 1], 
which represent all the transition rates of the system if 
the excluded volume effect between the motors is not 
considered. 
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where N indicates the microtubule site occupied by the 
leading motor, M is the total length of the membrane 
tube. Noting that N, M are not fixed numbers but vary 
with the process of tube growth, is the number of 
bound motor at microtubule site i, which is either or 1. 
rrij is the number of free motor at membrane site j . mo is 
the number of motor at membrane site 0, and P2N+3M+1 
represents the diffusion rate of the motor from membrane 
site (i.e., the reservior) to site 1. 

P2i-i corresponds to the unbinding rate of the mo- 
tor from microtubule site i, P2i corresponds to the for- 
ward stepping rate of the motor at microtubule site i. 
P2N+3j-2, P2N+3j-i Correspond respectively to the back- 
ward and the forward diffusion rate of a free motor at 
membrane site j. P2N+3j corresponds to the binding rate 
of a free motor at membrane site j to the microtubule. 

Considering the excluded volume effects, some of the 
above transitions are actually prohibited, so the corre- 
sponding transition rates P/ must be 0, or else P/ = P/. 
Especially, P2N+3M-1 = ^ because the free motors at 
membrane site M can not diffuse forward. 
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FIG. 2: (a) The threshold force Fm for membrane tube ex- 
traction as a function of the surface motor density po on the 
vesicle. The error bars represent standard deviation calcu- 
lated from 30 repeated simulations, (b) Velocity of membrane 
tube growth as a function of po when the retreating force is 
set as 27.5 pN. The error bars represent standard deviation 
calculated from 20 independent simulated membrane tubes. 
The simulation parameters are listed in Table |T] 



We define the global transition rate as 



2Ar-|-3M-|-l 



s= E Pi- 



(3) 



1=1 



In Gillespie algorithm, the time step At is not fixed and 
is a stochastic variable distributed exponentially with a 
characteristic time scale 1/S. At determines how long 
to wait to see a motor (no matter on which site) per- 
forming any one of the above mentioned transitions. 
To determine which transition will take place in the 
next time step, we generate a random number R dis- 
tributed uniformly in the range [0, S"]. The largest value 
of m e [1, 2iV + 3A/ + l] that fulfills the following inequal- 
ity, 



(4) 



1=1 
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means that the m-th transition numbered in Eq. [2] wih 
occur in the next time step. 

The simulation is performed as follows. The initial 
length of the membrane tube is set as M = 20 with each 
site occupied by a bound motor, and the actual transition 
rates P/ are calculated. Then the following steps are 
repeated. 

(1) Calculate the global transition rate S, and deter- 
mine the stochastic time step At. 

(2) Determine which transition to take place. In this 
step, one of the following four possible cases may be en- 
countered. 

a) A bound motor unbinds from the microtubule. So 
the number of free motors at the corresponding site of the 
membrane tube will increase by one. If this unbinding 
motor is a pulling motor, the load will be redistributed 
by other motors according to the equilibration of force, 
and the length of the membrane tube will also be changed 
too. 

b) A bound motor steps forward. Its tail domain will 
be dragged forward. If it is a pulling motor, the load will 
be redistributed among pulling motors, and the length 
of the membrane tube will be changed too. Sometimes 
an unpuUing motor can become a pulling motor after a 
forward step, as mentioned in Sec. |TT1 

c) A free motor diffuses from one membrane site to 
another. 

d) A free motor at membrane site i binds to the mi- 
crotubule, then m'^ = rrii — 1, Hi — 1. 

(3) After one of the above four cases performed, the 
array P/ is updated with the new configuration. 

(4) The time t is updated to t + At and the process is 
repeated. 

The repetition stops when it reaches the conditions 
that we set up, which will be explained in the following. 

At a constant motor density po, there exist a threshold 
retreat force Em above which the motors cannot extract 
a membrane tube from the vesicle. The Fm is determined 
as follows T7J. For a given value of the motor density po, 
we initially set the retreating force a very large value and 
repeat this process 200 times. If the force is too large, 
the membrane tube may retreat completely. If membrane 
tube retreat can be observed in all the 200 cases, we lower 
the force and repeat the process again till a threshold 
value Fjn is reached, at which at least one among the 200 
membrane tubes does not retreat. Typically, membrane 
tube retreat usually occurs when M < 50, and seldom 
occurs when M > 100. In our simulation, if the length 
of the membrane tube exceeds M = 250 (i.e. 2/xm), the 
membrane tube can grow persistently and no retreat can 
be observed, so we set M — 250 as the threshold value 
which means the forming of a membrane tube. 

To calculate the growth velocity of membrane tube, we 
let the membrane tube grows to a long enough length, 
6.4 /im. Since the system is not in steady state at the 
beginning, we calculate the average growth velocity when 
the length of the membrane tube is in the range from 0.4 
to 6.4 /im. 



IV. RESULTS AND DISCUSSION 

A. The growth of membrane tube 

In Ref . 0, , experiments showed that there exists 
a threshold value of motor density below which the mo- 
tors cannot extract membrane tubes from the vesicle at a 
given extraction force F. It means that at a constant mo- 
tor density pQ, the motors cannot pull a membrane tube 
if the retreat force is larger than a threshold force Fm- 
If Fm-po relation is known, we can predict the threshold 
density above which membrane tubes can be extracted 
from vesicles. In order to compare our results with Ref. 
prj . we use the same parameters as Ref. 17 1, which are 
listed in Table. HI And we also adopt the same method 
to find the threshold force as Ref. 17[ has done, as 
mentioned in Sec. IIIII 

In Fig. ma), we plot the simulated Fm-po relation. The 
force F„i increases with po from 14,pN (corresponding 
to po — 1-0 /ini"^) to 36 pN (corresponding to po = 
1000 pvcT^). In the experiments, the actual extraction 
force for membrane tube is found about 27.5 ± 2.5 pN 
Q. From Fig. HJa), we can know the corresponding 
threshold value of po to extract membrane tubes from 
vesicles is about 100 /xm~^, which agrees quite well with 
the experiment values about 100~200 /Ltm~^ H, Ell- 
in Table. HI all the parameters but the binding rate 
fcb, are adapted from experiments, fcb = 4.7 s~^ is a 
fitting value from Ref. [g, and its reliability would be 
questioned. In Fig|3l we plot the fcb-dependence of Fm 
at po = 200 pm"^. Fm increases fast with k^, when fcf, 
is smaller than 4.0 s^^, and then increase slowly to a 
saturated value about 36 pN. When > 2.0 s~^, P,„ will 
be larger than 27.5 pN, which is big enough to extract 
membrane tubes from a vesicle. 




FIG. 3: The threshold force Fm for membrane tube extraction 
as a function of the binding rate of motor fej, at the surface 
motor density po — 200 pm~^. The error bars represent stan- 
dard deviation calculated from 30 repeated simulations. 
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TABLE I: Parameters used in the simulations; reference sources are indicated. 



Parameter 


^spring 




D 


kb 


fcuO 




Value 


0.3 pN/nm 


3.0 pN 


1 pm^s"^ 


4.7 


0.42 


lOfcflT 


Ref. 


fl9] 




re, 17] 


[6,12 


fl7, 20] 


f6, 17] 



The growth velocity of membrane tubes as a function 
of po, at the given load of 27.5 pN, is shown in Fig. [2][b). 
The mean velocity is about 80 ~ 90 nm/s, consistent 
with the experiment value 90 ± 60 nm/s. The almost 
independence of the mean velocity on po is a prediction 
which can be tested by further experiments. 



B. The distribution of motors 

In Ref. [6], the motors are found not uniformly dis- 
tributed along the membrane tube. At the tip region, 
the density is bigger than that of the rest of the mem- 
brane tube. Our simulation results are in agreement with 
this observation. Fig. |3] shows that the mean density 
of motors increases from the vesicles to the tip of mem- 
brane tube. Because the motors in the tip region bear the 
load, their velocities are smaller than the motors in the 
rest region, so the motors will accumulate in the tip re- 
gion and generate enough force to extract the membrane 
tube. Especially at the very tip of the membrane tube, 
almost each microtubule site is occupied by a bound mo- 
tor. Moreover, there are free motors on the membrane 
tube, so the motor density in this region exceeds 1 per 
site. Since the motors in the tip region unbind from mi- 
crotubule more frequently than that of other region due 
to high load, the density of free motors in the tip re- 
gion is also higher than that of other region as shown in 
Fig. m These free motors will diffuse and rebind to the 
microtubule, and then step toward the tip region again, 
forming a cycle of motor flow. 

At the tip region of membrane tube, the number of 
pulling motors is also a stochastic variable which fluctu- 
ates during membrane tube growth. We now study the 
number of pulling motors and its probability distribution. 
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FIG. 5: (Color online) (a) The probability distributions of the 
number of pulling motors for two different motor densities, 
po ~ 200 pm~^ and po = 1000 pm~^ , respectively. The data 
are obtained from 10000 examples, (b) Average number of 
pulling motors function of time, which are averaged 

from 100 independent simulated membrane tubes. 
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FIG. 4: (Color online) The distribution of motors (i.e., the av- 
erage motor number per site) along the membrane tube. The 
data are averaged from 100 independent simulated membrane 
tubes at po = 200 pm~^. 



In Fig. [Sja) we plot the probability distributions of 
pulling motor number for two different motor densities, 
po = 200 pm~^ and po = 1000 /im~^ respectively. The 
data are obtained from 100 independent simulated mem- 
brane tubes, and 100 samples are picked out from each 
membrane tube at 0.5-second intervals during its growth, 
so there are 10000 samples for either po- Since succes- 
sive samples are separated sufficiently long time, they 
can be regarded as independent samples, so we can plot 
the histogram Fig. EJa) . During membrane tube growth, 
the average number of pulling motors is 8.4 and 8.5 for 
Po = 200 pm~^ and po — 1000 /zm~^ respectively. The 
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average number of pulling motors can be regarded as a 
time-independent constant, as shown in Fig. [5l[b). So the 
number of puling motors is almost independent of motor 
density and the length of membrane tube. On average, in 
the growth of membrane tube, each pulling motor bears 
3 ~ 4 pN, which smaller than the stall force of kinesin, 
about 7pN. 

V. SUMMARY AND CONCLUSION 

In this paper the mechanism of extracting membrane 
tubes by motors is studied. We propose a two-track- 
dumbbell model, by considering the excluded volume ef- 
fect of both the head domain and tail domain of motors, 
as well as the elasticity of the motor stalk. We conclude 
that (1) the membrane tube can be extracted by a group 
of motors along a single microtubule protofilament; (2) 
the growth velocity of the membrane tube is almost inde- 
pendent of the motor densities on the vesicle po! (3) the 
density of motors in the tip region is much higher than 
that in the root region. 

(2) is a prediction, (3) agree with experimental data 
0, and (1) is in contrary to the result of Ref. [13] which 
suggested that motors use several protofilaments simul- 
taneously to pull a single tube, without consideration of 
the structure and elasticity of the motor molecule. While 
these properties, especially the elasticity of motor stalks, 
can affect the multi-motor transport [141, it's also under- 
standable they play a role in extracting membrane tubes 



by a collect of motors. Our simulation results show that 
the number of pulling motors (which apply force against 
the retraction of membrane tube) at the tip region of 
membrane tube is about 8^9 and thus they can gen- 
erate enough force to extract a membrane tube from a 
vesicle and maintain its growth, while Ref. indicates 
the number of pulling motors is about 3 and thus three 
protofilaments are necessary to extract a membrane tube. 

Another difference between our model and Ref. [l7| . 
the choice of force- velocity {F-V) relation of a single mo- 
tor, may influence the simulation results too. In Ref. 
[l7| , the (F-V) relation of a single motor is simplified as 
linear, while we adopt the more reliable theoretical rela- 
tion from Ref. [ll| . In our previous study, we have shown 
that the property of multi-motor transport strongly de- 
pends on the F-V relation of a single motor [l3|- Hence, 
whether the membrane tube is pulled by one row or sev- 
eral rows of motors, and how the intrinsic properties of 
kinesin motors (the elasticity of kinesin stalk, the ex- 
cluded volume effect, and the characteristic F-V relation 
of a single motor) are involved in the extracting process, 
should be investigated by future experiments. 
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